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ABSTRACT 



Context. Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. In magnetic waveguides resonant absorption 
due to plasma inhomogeneity naturally transfers wave energy from large-scale motions to small-scale motions. In the cooler parts of 
the solar atmosphere as, e.g., the chromosphere, effects due to partial ionization may be relevant for wave dynamics and heating. 
Aims. We study resonant Alfven waves in partially ionized plasmas. 

Methods. We use the multifiuid equations in the cold plasma approximation. We investigate propagating resonant MHD waves in 
partially ionized flux tubes. We use approximate analytical theory based on normal modes in the thin tube and thin boundary approx- 
imations along with numerical eigenvalue computations. 

Results. We find that the jumps of the wave perturbations across the resonant layer are the same as in fully ionized plasmas. The 
damping length due to resonant absorption is inversely proportional to the frequency, while that due to ion-neutral collisions is in- 
versely proportional to the square of the frequency. For observed frequencies in the solar atmosphere, the amplitude of MHD kink 
waves is more efficiently damped by resonant absorption than by ion-neutral collisions. 

Conclusions. Most of the energy carried by chromospheric kink waves is converted into localized azimuthal Alfven waves that can 
deposit energy in the coronal medium. The dissipation of wave energy in the chromosphere due to ion-neutral collisions is only 
effective for high-frequency waves. The chromosphere acts as a filter for kink waves with periods shorter than 10 s. 

Key words. Sun: oscillations — Sun: atmosphere — Magnetohydrodynamic s (MHD) — Waves 



1. Introduction 

Recent observations have shown the ubiquitous presence of 
propagating magnetohydrodynamic (MHD) Alfvenic waves 
in the solar atmosphere. For example, Alfvenic transverse 
waves propagating in magnetic waveguides o f the solar 
corona were first ob served by iTomczvk et alj ('2007) and 
iTomczvk & Mc lntoshl (2009) using the Coron al Multichannel 
Polari meter (CoMP), and more recently by iMcIntosh et al. I 
(12017) using SDO/AIA. In chromospheric spicules, the presence 
of Doppler oscillation s is known f or more than 40 years (see the 
review bv lZaqarashvili & Erdelv i 2009). Recent observations of 
Alfvenic transver s e wav es in spicules have been rep o rted by, e.g, 
iDe Pontieu et alJ (l2007h IZa Qarashvi li"etan (l2007h. iKim et al. I 
( 20081) . iHe et al. I (l2009aibF . lOkamoto & De Po ntieu ("20111). In 
addition, propagating transverse waves in thin threads of so- 
lar prominences have been observed (e.g., iLin et aP l2007l 



l2009i) and Alfve nic waves in bright points have been reported 
( Jess et ai] l2009l) . The role and implications of the observed 
Alfvenic waves for the heating of the solar atmosphere have been 
discus s ed by, e.g., Erd elvi & Fedun (2007); De Pontieu et al] 
(l2007l) : IMcIntosh et al. I (l201 li) : ICargill & De Moortell (l20nl) . 

Based on MHD wave theory, a number of works have 
interpreted th e observed waves as p ropagating kink MHD 
waves (e.g., Erdelvi & Fedun 2007; Van Doorsselaere et al.' 
2008; Lin et al. 2009; Pa scoe e t al. 2010, 20 11; Terradas et aL 
20T0l IVerth et all l2010l 1201 ll: ISoler et aP IMlalblldh . Kink 
waves are transverse waves with mixed fast MHD and Alfvenic 
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prope rties (see, e.g., lEdwin & RobertsI 119831; iGoossens et al.l 
2009). In thin magnetic tubes kink waves are highly Alfvenic 
bec ause their domin ant restoring force is magnetic tension 
(Go ossens et al.ll2b09h . It has been shown that resonant absorp- 
tion, caused by plasma inhomogeneity in the direction trans- 
verse to the magnetic field, is a natural and efficient damp- 
ing mechanism for kink waves (see the recent reviews by 
IGoossens et al]|2006l; lGo"ossensll2008t iGoossens et alj|201 1*). In 
magnetic waveguides resonant absorption transfers wave energy 
from transverse kink motions to azimuthal motions localized in 
the inhomogeneous par t of the waveguide. This pr ocess has been 
studied numerically bv lPascoe et al.l (1201 Oll20 lib in the case of 
driven kink waves in coronal waveguides. 

Using analytical theo ry based on the thin tube and thin 
boundary approximations. iTerradas et alJ (1201 Ol hereafter TGV) 
obtained that the damping length due to resonant absorption 
is inversely proportional to the wave frequency. Therefore, it 
was predicted that high-frequency kink waves beco me damped 
in leng th scales shorter than low-frequency waves. IVerth et al] 
(1201 ol) showed that this result is consist ent with the CoMP 
observations of damped coronal waves ( Tomczvk et al J 120071 : 
Tomczvk & Mcintosh 2009). Subsequent investigations have ex- 
tended the original work by TGV by incorporating effec ts not 
considered in their paper. For example, 'Soler et al.l J2011b) took 
the presence of flow into account, and Soler et al. (201 Ic) stud- 
ied the influence of longitudinal density stratification. Both 
works concluded that the damping length remains inversely pro- 
portional to the frequency when flows and longitudinal stratifi- 
cation are present. 
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In TGV and subsequent works cited above, the plasma 
is assumed fully ionized. However, the plasma in the cooler 
parts of the solar atmosphere is only partially ionized as, e.g., 
in the chromosphere or in prominences. This fact raises the 
relevant question on whether the damping length remains in- 
versely proportional to the frequency when the plasma is par- 
tially ionized or, on the contrary, this dependence is modi- 
fied by the effect of ion-neutral collisions. The effect of par- 
tial ionization on the damping of Alfven waves has been in- 
vestigated in a large number of papers in differen t contexts. 



Some examples are the works by, e.g., iHaeren del 



(19921); 

PecseH & Engvold (2000); 
(2007); Soler et al.. 



De Pontieu & Haerendel (1998); 

De Pontieu et al, (2001) ; Fortez a et al. , . . , 

(2009a); Car boneU et all (l2010l) ; Izaqarashvili etanl201 lalbl) . 
However, these works only focus on the role of partial ioniza- 
tion for the damping and do not consider resonant absorption, 
which is a basic and unavoidable phenomenon when plasma 
and/or magnetic inhomogeneity is present. 

To our knowledge, the first atte mpt to study resonan t waves 
in partially ionized plasmas was bv Soler et all (l2009b^. These 
auth ors used the single-fluid approximation (see, e.g., Braginskii 
1965 ) to investigate standing resonant kink waves in a model of 
a partially ionized prominen ce thread. Subseque ntly, this first in- 
vestigation was extended in ISoler et al] (1201 lai hereafter SOB) 
to the case of propagating waves. The use of the single-fluid ap- 
proximation, as in Soler et al. (2009b) and SOB, seems reason- 
able as the expected values of the collision frequencies in the so- 
lar atmosphere are much larger than the observed wave frequen- 
cies. However, the single-fluid approximation misses important 
effects when small length scales and/or high frequencies are in- 
volved. In such cases, th e multifluid description is a more suit- 
able approach (see, e.g., IZaqarashvili et al.ll201 lallbh . For reso- 
nant waves, perturbations develop very small length scales in the 
vicini ty of the resonance p ositio n (see, e.g., [Ti rrv & Goossens 
1996; Ruderman & Wright" 1999t IVasquez(l2005l: iTerradas et al.1 
2006), where multifluid effects may play a relevant role. 

Here we perform a general description of propagating reso- 
nant Alfvenic waves in partially ionized plasmas using the mul- 
tifluid treatment. Our work extends the investigation of SOB by 
considering arbitrary values of the collision frequencies and by 
performing a more in-depth analysis of the resonant process. 
Section |2] contains a description of the equilibrium configura- 
tion and the basic equations. In Section [3] the behavior of wave 
perturbations around the Alfven resonance position are investi- 
gated and jump relations for the perturbations are derived. Later, 
in Section |4] we perform an application to resonant kink waves 
in straight tubes and obtain expressions for the damping lengths 
due to resonant absorption and ion-neutral collisions. The ap- 
proximate analytical theory is complemented with fully numer- 
ical eigenvalue computations. In Section |5] the implications of 
our results for the particular case of chromospheric waves is dis- 
cussed. Finally, we give our main conclusions in Section|6] 

2. Equilibrium and basic equations 

We consider a partially ionized plasma composed of ions, elec- 
trons, and neutrals. We use cylindrical coordinates, namely r, ip, 
and z for the radial, azimuthal, and longitudinal coordinates. The 
medium is permeated by an equilibrium magnetic field, B, of the 
form 



B = B^e^ + B-e 



(1) 



In general, both azimuthal, B^, and longitudinal, B., components 
are functions r. The equations governing the dynamics of a mag- 



neti zed multifluid pl asm a are discussed i n classical works as, 
e.g.. Cowling 71956) and Braginskii' (1965), or more recently in 
Birk et al. ( 1996) and Pinto et al. (2008). Recent investigations 
which make extensive use of the m ultifluid descript ion in the 
context of MHD waves are Zaqarashvili et al.l (1201 lallbl) . We re- 
fer the reader to these two papers for details about the derivation 
of the equations. In brief, a partially ionized multifluid plasma is 
governed by the equations of the different species, which contain 
terms that couple the various fluids by means of collisions. 

Here we study linear perturbations superimposed on the 
equilibrium state. Due to the very small momentum of electrons, 
we neglect the collisions of electrons with ions and neutrals. In 
addition, we consider a magnetically dominated plasma and ne- 
glect the gas pressure of ions compared to the Lorentz force, and 
the gas pressure of neutrals compared to the collisional friction 
with ions. This simplification neglects the plasma displacement 
along the magnetic field direction and removes the slow or cusp 
continuum from the equilibrium. Thus, in the present work we 
focus on resonance absorption in the Alfven continuum only. 
The inclusion of gas pressure is a subject for future works. Under 
these conditions, the basic equations of our investigation are 



Pn 



5v, 

' dt 

dt 
db 

dt 



1 



(V X b) X B - PiPnJin (Vi - Vn) , 



-PiPnTin (Vn - Vi) , 

V X (vi X B) , 



(2) 
(3) 
(4) 



where Vi and ¥„ are the velocities of ions and neutrals, respec- 
tively, b is the magnetic field perturbation, pi and p„ are the den- 
sities of the ion and neutral fluids, respectively, yU is the magnetic 
permittivity, and is the ion-neutral collision rate coefficient 
per unit mass. Instead of using in the remaining of this paper 
we use the ion-neutral collision frequency, vin, which has a more 
obvious physical meaning. The ion-neutral coUision frequency 
is defined as 



Vin = PiTin 



(5) 



Equations (|2]i and ([3]) are the linearized momentum equations 
of ions and neutrals, respectively. Equation (HJi is the linearized 
induction equation, in which we have omitted all nonideal, dif- 
fusion terms as, e.g., magnetic resistivity. We do so because the 
goal of the present paper is to assess the particular role of ion- 
neutral collisions on the resonance absorption process. The ef- 
fects of other dissipative mechanisms as, e.g, resistivity or vis- 
cosity, have been extens ively studied in the e xisting literature 
(see the recent review bv lGoossens et alj|201 lb 

In this work the equilibrium quantities are functions of r 
alone, so that the equilibrium is uniform in both azimuthal and 
longitudinal directions. Hence we can write the perturbed quan- 
tities proportional to exp (invp + ik~z), where m and k, are the az- 
imuthal and longitudinal wavenumbers, respectively. We express 
the temporal dependence of the perturbations as exp(-;wf), 
with (jj the frequency. Since there are no equilibrium flows, the 
Lagrangian displacement of ions, ^i, and neutrals, ^n, are 



to 



(6) 



From Equation ^ it is straightforward to derive the relation 
between ^n and ^i, namely 



U) + iVin 



(7) 
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Equation (|2| informs us that and are related to each other. 
The factor of proportionally is a complex function of the wave 
frequency, uj, and the ion-neutral collision frequency, Vin. Hence, 
there is a phase difference between the motions of ions and neu- 
trals. Physically, this phase difference can be interpreted as the 
"time delay" between the motions of ions and the corresponding 
reaction of neutrals. Thus, for vin s> w, = and the phase 
difference is zero, i.e., ions and neutrals behave as a single fluid. 
On the contrary, for Vin «c u ions and neutrals become decou- 
pled from each other, i.e., the coUisionless case. In the remaining 
of this paper, we call the limits Vi,, » w and Vin <K w the single- 
fluid limit and the coUisionless limit, respectively. 

Although we perform a general analysis for arbitrary vin and 
a), we must note that vin » w is the realistic situation according 
to the observed frequencies of MHD waves in the solar atmo- 
sphe re and the expected value s of the colhsion frequency (see , 
e.g., iDe Pontieu & Haerendell 119981: iDe Pontieu etal] l2001h . 
Thus, the case Vin » w will receive special attention. 

3. General expressions for resonant Alfven waves 

In this section we study the properties of resonant MHD waves 
in the Alfven continuum. We follow the notation used in, e.g., 
ISakurai et al] (IT991): Goo ssens et al.1 ([1992, 1995). We com- 
bine Equations (l2]l-(|4]i to arrive at two coupled equations for 
the radial component of the Lagrangian displacement of ions, 
^ir = /Vi ■ Crjoj, and the Eulerian perturbation of the total pres- 
sure, f = B ■ hj^i, namely 



dr 

2)^ = Cs^ir - CiP, 
dr 

with 
Ci 



IB, 



0) 



2 m 



Id r 



0?- 



V 



m 



C3 = © 



R 

. - ^ 2 2 f 

+ ——(^ - -piVAi^A — ' 

yi ju 



(8) 
(9) 

(10) 
(11) 
(12) 

(13) 



where v^j - ^ is the square of the Alfven velocity of the ions, 
with - B\ + Z??. In addition we have defined 



3A. Modified Alfven continuum 

Equations (|8]l and (|9]i are singular when D = 0. As the equilib- 
rium quantities are functions of r, the position of the singularity 
is mobile and depends on the frequency, co. The whole set of 
frequencies satisfying D - at some r form a continuum of 
frequ encies called the Alfven continuum (see, e.g..l Appert et al.l 
Il974l) . In the ideal, fully ionized case the condition D — is 
satisfied where the frequency, oj, matches the local Alfven fre- 
quency, LOa- 

We investigate how the Alfven continuum is affected by ion- 
neutral collisions. The equation 2D = 0, which D defined in 
Equation ( fTOl i. can be rewritten as a third-order polynomial in 
w, namely 



o)^ + i( l + a) ViaLi/ - coicj ■ 



0. 



(18) 



Equation (fTST l describes the Alfven continuum modified by ion- 
neutral collisions. For vin = we recover the ideal Alfven con- 
tinuum in a fully ionized plasma, i.e., oj - +a)A- Note that 
Equation (fTST l and so the continuum frequencies are independent 
of m. 

To shed light on the effect of ion-neutral collisions, it is con- 
venient to rewrite Equation (fTSl l in the following form 



2 2 



1 + /i 



1 4-/(1+ a) 2^ 



= 0. 



(19) 



This enables us to easily evaluate the continuum frequencies at 
the limits vin » o) and Vin <k lo, namely 



+(L)A (1 + ») , for Vin » W, 
+CiJa, for Vin ^ OJ- 



(20) 



In the limit Vin <s; w we obtain again the continuum frequencies 
of the ideal, coUisionless case. In the limit Vin » w, the contin- 
uum frequencies are modified by the plasma ionization degree 
but the frequencies remain real, so that the position of the sin- 
gularity of Equations (O and (|9]l is shifted with respect to the 
coUisionless case. 

The limits Vin » co and Vin <k w are extreme cases. Now, 
we explore the continuum frequencies for arbitrary Vin to have 
a complete picture of the role of ion-neutral collisions. For ar- 
bitrary Vin we find an approximate solution to Equation (fTST l by 
assuming that the effect of collisions on the continuum is to pro- 
duce a weak damping of the continuum modes. This means that 
the continuum frequencies are assumed complex and that their 
real parts are taken the same as in the coUisionless case. Thus 
we write w * +cl>a + ij, with y <K aj\, and put this expression in 
Equation (fTSl l. After neglecting terms with 0(y^^ we obtain the 
modified Alfven continuum frequencies, namely 



/b 



-B^ + k.B,, 



+ / (I + a) Vif 



Pn 

Pi ' 



(14) 

(15) 

(16) 
(17) 



In these expressions, wa is the local Alfven frequency for the 
ions, CJ is the effective or modified frequency due to colhsions, 
and a indicates the plasma ionization degree. 



1 

+Ci>A - l^aVin. 



(21) 



Equation (|2TJ shows that the continuum frequencies in a par- 
tially ionized plasma are complex. Based on this result we antic- 
ipate that waves driven with a real frequency a> will not produce 
true singularities in Equations ([8]l and ^ since the Alfven con- 
tinuum frequencies are in the complex plane. This means that the 
singularity is removed in a partially ionized plasma even in the 
absence of additional dissipative mechanisms as, e.g., resistivity 
or viscosity. 

As the two solutions given in Equation ( |2TI ) are complex, 
the third solution to Equation ( fTSb must be purely imaginary. 
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So, we write oj x iy and neglect terms with O (y^) to obtain an 
expression for the third solution, namely 



U K —IV., 



(22) 



which corresponds to a purely damped, nonpropagating mode 
whose damping rate is given by the collision frequency. The 
presence of this third solution is not relevant for our subsequent 
analysis. 

3.2. Behavior of perturbations around tiie resonance position 

Here our purpose is to assess the behavior of the perturbed quan- 
tities around the position of the ideal Alfven resonance. We de- 
note by ta the ideal Alfven resonance position. The value of is 
obtained from the ideal resonance condition, i.e., = w^(rA), 
where a> is assumed to be real. Around r = the behavior of 
the perturbations is dominated by the resonance. At r = ta ion- 
neutral collisions become relevant to remove the singularity of 
the coefficient D in Equations ^ and while their effect on 
the rest of coefficients is minor. This means that around r - rp, 
we can take the ideal expressions of coefficients C\, Ci, and C3, 
and only keep the terms related to collisions in the expression of 
coefficient D. 

We define the new variable s - r - rp^. The ideal coefficients 
Ci, C2, and C3 at = become the following constant terms. 



C 



2 B^B~ /b^b 



4 BlBl fl 



C3 



(23) 
(24) 
(25) 



with gB = f^B^ - k,B^. All quantities in Equations (l23Tl-(l24li 
have to be evaluated at i = 0. Then we can rewrite Equations dSjl 
and (|9]l as 



D 



ds 



di 
dCA 



D 



^Ca, 

2 B^B, /b 
fl Ta fipi 

0, 



Ca 



ds 

where Ca is defined as 



Ca^SbP- — 



2 B^B, 



/" '"A 



(26) 
(27) 
(28) 

(29) 



The quantity Ca is approximately conserved across the res- 
onance position. This is the same cons ervation law a s that 
in ideal and dissipative MHD (e.g .. Sakurai et alJ Il99lt 



m Ideal and dissipative MhUJ (e-g-- laaJcurai et alJ 11^^ it 
iGoossens et al.|[T995t iTirry & Goossensll 19961) . The Lagrangian 
displacement perpendicular to the magnetic field lines, - 
(^i^Bj - ^izB^) /B, is related to Ca as 



B ^ 
i— Ca. 



(30) 



Next we approximate the coefficient D by its first-order 
Taylor polynomial around s - 0. The linear expansion of D is 



approximately valid in the interval -^a < * < ■sa, with sp^ satis- 
fying the relation ^a |(w^)7(<j^a) "I' where the prime denotes 
radial derivative. The resulting expression is 



D ~ PiviiAA (A H- s) , 



with 



ds 



(31) 



(32) 



and A - Ar + /Aj, where Ar and Aj are the real and imaginary 
parts of A, respectively, given by 



Ar 



Aa o)^ + 



A, 



Aa aP- + 



(33) 



As before all quantities in the previous expressions have to be 
evaluated at i = 0. We find a constant complex term. A, in the 
expansion of D (Equation (ISTTi). The presence of this term re- 
moves the singularity of the solutions, so that D Q for s - 0. 
We use Equation dJTT i to rewrite Equations (l26T l and (l27T i as 



d^ir _ gB Ca 1 
ds fipi^v^. Aa ^ + s' 
dP IB^B, /b Ca 1 



ds n rA iip{^v\^ Aa A + s' 



(34) 
(35) 



For our subsequent analysis it is convenient to use the scaled 
variable t defined as 



s + K\ 



R 



Ai 

Then we integrate Equations ( l34l ) and ( l35T l and obtain 

/^prv^j Aa 

2 Ba;B, fp> Ca 

P ^ - ^ ' -^0(t) + Cp, 

fi rA ijpih'j^.AA 



(36) 

(37) 
(38) 



where C^ and Cp are constants of integration. The ^(t) function 
is defined as 

g(T) = ln(T + i). (39) 
We separate the real and imaginary parts of ^(t), namely 

(40) 



^ (&(t)) = In Vr2 + 1, 



n ] arctanr ' -htt, for r < 0, 

^(^^^^^^^arctanr-', for r > 0. 



(41) 



Note that we need to add n to the imaginary part of Q(t) for 
T < in order to make the solution continuous and analytic at 
T = 0. This is not a problem since we can always incorporate a 
constant of integration. 

From Equation (|30] | we obtain that the perpendicular dis- 
placement behaves as 



Ca 



1 



-T{t\ 



PiB AaAi 
with T{t) defined as 
/ 



T{t) = 



T + I 



(42) 



(43) 
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(a) 1.0 



0.5 - 



0.0 



-0.5 - 



-1.0 





-10 



Fig. 1. Real (solid line) and imaginary (dashed line) of functions 
a) T{t) and b) ^(t). 



We separate the real and imaginary parts of 'F(t), namely 



3(^r(T)) 



t2 + 1 



(44) 



the 



Functions Tir) and Q{t) play here the role of 
universal functions found in a number of previous inves- 
tigations of resonant waves in both stationary and non 
stationary 
(see, e.g.. 



st ates and for differen t dissipative pr o cesses 
Mok & Einaudi' '1985; 'Goossens et al.' '\99^, 



Ruderman et al 
1991 



1995; 



Tirrv & Goossens 1996; Erdelvi et al. 
IWright & Allai| Il996t iRuderman & Wrightl " Il999t 
32). A comprehensive review on the 
T{t) and Qij) 



Vanlommel et al. 200: 



A comprehensive 
properties of the universal 



importance and 

functions can be found in Goos sens et alj (1201 ll) . In particular, 
ou r TjT) and Q{t) f unctio ns coincide with the functions found 
by IWright & AllanI (1 19961) for magnetospheric Alfven waves 
damped by Pedersen conductivity (see their Equations (23) 
and (24)), and are also equivalent to the functions described by 
IVanlommel et al. (2002) for non-stationary quasi-modes in the 
cusp continuum (see their Equation (16)). Figure [1] displays the 
real and imaginary parts of functions T{t) and 0{t). 

Ion-neutral collisions generate a collisional layer around the 
location of the resonance, i.e., t = 0, in a similar way as re- 
sistivity and/or viscosity generate a dissipative layer (see, e.g . , 
iHollweg & Yangll988HPoedts et alJl99"0al:ISakurai et al.lll99lh . 



The thickness of the collisional layer, 6, is proportional to Aj, 
namely 



6 ~ lAil = 



aajvir 



IAaI a;2 + vl 



(45) 



It can be seen that 5 — > in both limits Vin » oj and Vin <k 
o), whereas 6 takes its maximum value for vin - (d. This result 
points out that the collisional layer is extremely thin for realistic 
collision frequencies, i.e„ vin » w. 

3.3. Jump Conditions 

Here we use the expressions derived in the last section to deter- 
mine the jump of the perturbations across the collisional layer. 
In the scaled variable r and assuming a> > Q, the limits s ±s\ 
are equivalent to the limits t — > ±sign(AA)co. We denote by [X] 
the jump of the quantity X across the collisional layer, which we 
define as 



[X] = lim X{t) - lim X(t). 

T^sign(AA)c« T^-sign(AA)M 

Hence, the jumps of and P are 

Aa 

[P] = - 2 2 

IJ. rA m Aa 



(46) 

(47) 
(48) 



To calculate the jumps of ^i, and P we need the jump of 0, 
namely 



[g] = -/sign(AA)7r. 
Then it is straightforward to obtain 



[P] = 



8b Ca 

—171 

//p,2y2 . |Aa| ' 

. 2 B^Bz /b Ca 

— 171 

^l rA m^^ii IAaI 



(49) 

(50) 
(51) 



We find the remarkable result that these jumps are the same 
as obtained f or ideal and dissipative MHD (e.g., ISakurai et all 
1991; Goosse ns"etal] [T995; Ruderman et al] |1995l) . This means 
that ion-neutral collisions do not modify the jumps of the plasma 
perturbations across the resonant layer. In the particular case of 
a straight magnetic field, i.e., - 0, the conserved quantity 
across the resonant layer is proportional to the total pressure per- 
turbation, and the jump conditions for and P (Equations (l50l l 
and ( fSTl i) become 



np-/r\ 

T J 

Pi |AaI 



[P] = 0. 



(52) 



Finally, the asymptotic expansion of function T for t — > ±oo 
gives as the asymptotic behavior of when we move away 
from the resonance position as 



Ca 



1 



PiB AaAi t' 



(53) 



which again corresponds with the asymptotic dependence t"' 
found in both ideal and dissipative MHD. 

In this Section we have studied the case of propagating 
waves, i.e., real w and complex k,. However, most of the analysis 
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is equivalent for the case of standing quasi-modes, i.e., complex 
(jL> and real k,. In particular, the jump conditions (Equations (ISOl l 
and ( BTI )) are exactly the same for both propa gating and stand- 
ing modes. We refer the read er to, e.g.. lTirrv & Goossensl(ll996l) 
and IVanlommel et al.l (12002 ^ for details about the derivation of 
the jump conditions for standing quasi-modes. 



4. Kink waves in straight tubes 

Here we perform an application of the theory of Section|3]to the 
case of propagating resonant kink waves in straight flux tubes. 

We consider a straight magnetic cylinder of radius R with a 
constant vertical magnetic field, i.e., B^p - Q and B~ - constant. 
The tube is homogeneous in the longitudinal direction. In the 
radial direction, the tube is composed of an internal region with 
constant ion density pi i , a nonuniform transitional layer of thick- 
ness /, and an external region with constant ion density pi2- In the 
nonuniform layer the density changes continuously from the in- 
ternal to the external values in the interval R - 1/2 < r < R + 1/2. 
Hereafter, indices '1' and '2' refer to the internal and external 
regions, respectively. For the sake of simplicity, the parameter 
a is taken constant everywhere so that the neutral density, p„, 
follows the radial dependence of the ion density. The ion-neutral 
collision frequency, Vin, is also a constant to simplify matters. 

We restrict ourselves to the case of thin tubes, i.e., the wave- 
length is much longer than the radius of the tube. For propagat- 
ing waves the thin tube (TT) approximation is equivalent to the 
low-frequency approximation, i.e., (jjR/\\ «c 1, with \\ the kink 
velocity of ions defined as 



Pil +Pi2 



To check whether this approximation is realistic, let us 
consider the properties of the obs erved waves propagat- 
ing in chromospheric spicules (e.g., iDe Pontieu et alJ l2007t 
[Okamoto & De Pontieu 2011). Taking R = 200 km for flie 
spicule radius and Vk - 270 km s~' for the average d phase ve- 
locity as estimated bvl Okamoto &De Pontieu' ('2011'), we obtain 
toR/vi^ ^0.1 for a period of 45 s (Okamoto & De Pontieu 201 1) 
and cjR/vi^ ^ 0.02 for a period of 3 min (De Pontieu et al.ll2007ft . 
Therefore, the use of the TT approximation is justified. 

Addit ionally, we use the Thin Boundary (TB) approximation 
(see, e.g.. lHollweg & Yandl 19881) and assume l/R«: 1. The TB 
approximation enables us to use the jump conditions given in 
Equation ( |52] ) as boundary conditions for the wave perturbations 
at the tube bound ary. Detailed explanati ons of this method ca n 
be found in, e.g., iGoossens et al ] (l2006h andiG oossen 1 (I2008h . 
Doing so, we arrive at the dispersion relation for resonant kink 
(m = 1) and fluting (m > 2) waves in the TT and TB approxima- 
tions. For simplicity we omit the interm ediate steps which can 
be found in, e.g. . .Goossens et al.l(ll992l) . The dispersion relation 



IS 



m/R 



(55) 



In Equation ( ISST l we approximated ~ R. The effect 
of ion-neutral collisions is enclosed in the definition of 
CO (Equation (fTSIl). For fixed and real oj the solution of 
Equation ( l55b is a complex longitudinal wavenumber, - 
k,R + ik,i, where r and k^i are the real and imaginary parts of 
k^, respectively. In the present normal mode analysis, which rep- 
resents the stationary state of wave propagation, the amplitude 



of the propagating wave is proportional to exp (-z/Ld), with Ld 
the damping length defined as 



1 

Ln = — . 



(56) 



We investigate kink waves and set m - 1 . In the absence of 
resonant absorption the right-hand side of Equation ( |55] ) is zero. 
In such a case, we can obtain an approximate expression for k^^ 
by assuming k^^ <sc k^^, namely 



k^ 



0^ (A? + {l +a)vl 
•1 up- + y?„ 



(57) 



ur (1 + a) /v 



^, while 



iiP/vi- Note that the 



In the single-fluid limit, Vin » w and kP:^ 
in the fully ionized case a = and kr^ ~ ^ , 
expressions for k?:^ in the single-fluid and fully ionized cases dif- 
fer by a factor (1 + a). SOB studied propagating kink waves in 
the single-fluid approximation and found the same expression 
for fc^j, for both partially ionized and fully ionized plasmas (see 
their Equations (11) and (15)). The reason is that SOB used a 
slightly different definition of the kink velocity, Vk. Here, the 
kink velocity is defined using the ion density only, while SOB 
defined Vk using the total (ion H- neutral) density. In SOB the to- 
tal density is fixed and the effect of the ionization degree is to 
change the relative ion and neutral densities. On the contrary, 
here the ion density is a fixed parameter while the amount of 
neutrals depends on the ionization degree. 



(54) 4. J . Approximate expression for tlie damping lengtli 



Here we seek an approximate expression for Lu- First we evalu- 
ate the factor piiR) |AaIs in Equation ( l55l l by using the resonant 
condition aP - oj^ (R), namely 



Pi(R) \Aa\r = 



dpi 



dr 



(58) 



where |dpi/dr|^ is the radial derivative of the density profile at 
r - R. Next we write k, - k,j^+ik,i and assume weak damping, so 
we neglect terms with O (k^XAn expression for the ratio k,i/k,^ 
is obtained from Equation ( l55l l. namely 



\2 
■Pi2) 



1 



k-i 1 awvin 

k,R 2 oj^ + (I + a) v^^ SR Pii+Pi2 |dpi/dr|s" 
Now we express |dpi/dr|g as 



(59) 



dpi 
dr 



Pil -Pi2 

4 I ' 



(60) 



with F a numerical factor that depends on the form of the 
density profile. For example, F = A/n^ for a linear profile 
(Goossens et al. 2002 ) an d F - 2/n for a sinusoidal pro- 
file (Ru derman & RobertsI 12002). We assume that ^ r is ap- 
proximately the same as in case without resonant damping 
(Equation ( l57b ) and work on Equation ( |59] | to find the expres- 
sion for Ld as. 



1 



(61) 
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with Ld,ra and Ld.in the damping lengths due to resonant ab- 
sorption and ion-neutral collisions, respectively, given by 



Ld,ra = 2;r^Vk- 



1/2 



I (-Icj 



w2 +a)vl) 

1/2, , ,n1/2 



(62) 
(63) 



with f = Pii/Pi2 the ion density contrast. Importantly, we find 
that both the collision frequency, Vin, and the ionization de- 
gree, a, are present in the expression of the damping length 
by resonant absorption (Equation (l62b). In the single-fluid limit 
(vin » a>) Equations ( |62] ) and ( |63] | become 



-^D.RA 



■^D.IN 



2vk 



(l+«) 



1/2, 



1 



a 



(64) 
(65) 



The damping length due to resonant absorption in the single- 
fluid limit is inversely proportional to the frequency as in TGV. 
Indeed, in the fully ionized case {a - 0) Equations (|62] | and 
(l64l l become Equation (22) of TGV. For a partially ionized 
plasma Ld.ra also depends on the ionization degree, a. The de- 
pendence on the ionization degree was not discussed by SOB. 
SOB obtained an expression for the ratio of the damping length 
to the wavelength (see their Equation (20)) and did not ex- 
plicitly write the expression for the damping length. It turns 
out that the factor containing the ionization degree cancels out 
when the ratio of the damping length to the wavelength is com- 
puted, and so the dependence of Ld.ra on the ionization de- 
gree was not noticed by SOB. The damping length by ion- 
neutral collisions in the single-fluid limit (Equation ( |65]) ) is in- 
versely proportional to a P' . This result is co n sistent with the ex- 
pressions found by, e.g , ^Haere ndell (Il992h : iPecseli & Engvoldl 
( 120001) : iDe Pontieu & H aerendel d 19981) . and SOB among oth- 
ers. 

Figure |2] shows Ld/^ versus Vin/w for a particular set of pa- 
rameters given in the caption of the Figure. The total damping 
length (solid line) has a minimum for Vi,, ~ w. The reason for this 
behavior is that damping by ion-neutral collisions (dotted line) 
becomes more relevant than resonant damping (dashed line) for 
Vin ~ In both the collisionless (vin <s: o)) and single-fluid 
(vin » oj) limits, the total damping length is well approximated 
by the damping length due to resonant absorption. In particular, 
for Vin ^ the result tends to that predicted by Equation ( l64l l. 

Now we fix the collision frequency to VmRI\\ - 100 and 
compute LqIR as a function of (l)R/i\ (see Figure |3]l. The re- 
maining parameters are the same as in Figure |2] Although the 
observed wave frequencies in the solar atmosphere correspond 
to wTJ/vk «; 1, it is interesting to perform a general study 
for high values of the frequency. Our results point out that the 
damping of the kink wave is governed by diff'erent mechanisms 
depending on the value of wiJ/vk. We find that the damping 
length for ljR/v]^ « 1 is dominated by resonant absorption. Ion- 
neutral collisions start to become important when w/J/vk ~ 10. 
Finally, collisions are the dominant damping mechanism for 
large a)R/v\^. We must note that the analytical solution in the 
TT approximation may not provide accurate results for values 
of (L)R/i\ departing from the limit <jjRI\\ <«; 1. For this reason, 
numerical eigenvalue computations which overcome the limita- 
tions of the analytical approximations are performed in the next 
Subsection. 
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Fig. 2. LqIR vs. Vin/(jj for the kink wave in a magnetic flux tube 
with II R = 0.2, a = 0.5, f = 10, F = 2/7r, and w/?/vk =0.1. The 
solid line is the total damping length in the TT and TB approx- 
imations computed from Equation ( |6TI ). The dashed line is the 
damping length due to resonant absorption (Equation ( l62b ) only, 
and the dotted line is the damping length due to ion-neutral col- 
lisions (Equation (163)) ) only. Symbols O correspond to the full 
numerical eigenvalue results. 




Fig. 3. Ly^IR vs. uRjv^ with Vin^/vk - 100. The remaining pa- 
rameters and the meaning of the line styles are the same as in 
Figure |2] 



4.2. Numerical results 

Here we compare the approximate analytical result obtained by 
solving the dispersion relation in the TT and TB approxima- 
tions (Equation dHTT l) with the damping length obtained by solv- 
ing the full eigenvalue problem numerically. The numerical so- 
lution is not limited by the TT and TB approxim ations. The 
numerical code is similar to that used in TGV and ISoler et al.l 
(2011b). The reader is refereed to these papers for details of 
the numerical scheme. In short. Equations (|2])-(|4l) are integrated 
in the radial direction assuming a time dependence of the form 
exp {-iujt) and a spatial dependence in the tp and z directions as 
exp (itp + ik,z). The code solves the eigenvalue problem for the 
temporal damping of standing waves, i.e., complex u provided 
a fixed and real k^. To study spatial damping we need to convert 
the results from complex cj and real k, to real lo and complex k,. 
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The conversion is done following the method explained in TGV 
(see their Equation (40)). 

First, we consider the same parameters as in Figure |2] and 
compute the numerically determined Lx^jR versus Vin/w. To 
compare with the analytical approximation, the numerical result 
is plotted using symbols O in Figure |2] A very good agreement 
between approximate and numerical results is found. This means 
that for (jjRlv\^ <K 1 the approximate analytical theory provides 
accurate results. 

Next we numerically explore the effect of increasing the 
wave frequency. We take the same parameters as in Figure [3] 
Again, we use symbols O to represent the eigenvalue result 
in Figure [3] We find that in the eigenvalue computations the 
transition between the regime dominated by resonant damp- 
ing and that dominated by collisional damping occurs around 
(jjRI\\ ~ 1, while in the analytical approximation collisions start 
to become important for ojR/vi^ ~ 10. This discrepancy is an ef- 
fect of the TT approximation. For realistic values of the wave 
frequency (w/?/vk ~ 10"^ - 10"') both numerical and analytic 
results are in excellent agreement. 



5. Chromospheric kink waves 

The results of Section|4]have direct implications for MHD waves 
propagating in partially ionized plasmas of the solar atmosphere. 
For kink waves studied in this paper, both resonant absorption 
and ion-neutral collisions decrease the amplitude of the waves. 
However, the two processes represent very different physical 
mechanisms. 

On the one hand, resonant absorption is an ideal process that 
transfers wave energy from global kink motions to localized az- 
imuthal motions within the transversely inhomogeneous part of 
the flux tube. These azimuthal motions keep propagating along 
magnetic field lin es (se e the numerical simulations by, e.g., 
[Pascoe et al, 20ia l201lh . A detailed investigation of the energy 
transfer in the case of standing waves was done in lArregui et alJ 
(1201 1) by analyzing the Poynting flux in a two-dimensional con- 
figuration. However, resonant absorption itself does not dissipate 
wave energy in the plasma. The energy fed into the inhomoge- 
neous layer will b e dissipated by anot he r mech anism later (see, 
e.g., the results of iPoedts et af]ll989allbl IT990allblQ in resistive 
MHD). Hence, the damping length due to resonant absorption, 
Ld,ra, represents the length scale for the kink motions to be con- 
verted into azimuthal motions. 

On the other hand, ion-neutral collisions is a true dissipative 
process which deposits wave energy in situ and so it contributes 
to plasma heating. Hence, the damping length due to ion-neutral 
collisions, Ld.in, represents the length scale for the kink wave 
energy to be dissipated by ion-neutral collisions. By comparing 
the values of Ld.ra and Ld,in we can estimate the fraction of en- 
ergy converted to Alfvenic, azimuthal motions and the fraction 
of energy dissipated by collisions. 

Let us apply this theory to kink waves propagating along 
chromospheric waveguides (spicules). We assume that the driver 
of the waves is located at the photosperic level and the waves 
propagate through the chromosphere to the corona. We take 
the variation of physical parameters with height (e.g., density, 
temperature, ionizati on degree, etc.) from the VALC model 
dVernazza et al.|[T98l]) . For the chromospheric m agnetic field we 
consider the model used bv lLeake & ArbeJ (12006). Then we use 
Equations (|64] | and ( |65] ) to compute the values of Lq.ra and 
i'D.iN- As the physical parameters change along the spicule, both 
Ld.ra and Ld.in are functions of height in the chromosphere. We 
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Fig. 4. (a) Ld.ra (solid line) and Ld.in (dashed line) vs. height 
above the photosphere for a wave period of 45 s. (b) Averaged 
Ld.ra (solid line) and Ld.in (dashed line) in the chromosphere 
vs. wave period. The horizontal dotted line represents the height 
of the chromosphere above the photosphere. 

plot in Figure |4fa) the values of the damping lengths as func- 
tions of height for a wave period of 45 s (Okamoto & De Pontic^ 
1201 1). The damping length due to resonant absorption increases 
with height. Th is is an effect of th e increase of the kink velocity, 
v'k, with height (ISoler etalJ2011cl) . At low heights, Ld.ra is com- 
parable to the thickness of the whole chromosphere, meaning 
that a large fraction of wave energy is in the form of azimuthal 
motions when the wave reaches the coronal level. On the con- 
trary, Ld.in is several orders of magnitude longer As Ld,in de- 
creases with height, the eff'ect of collisions is more important in 
the upper chromosphere. 

Next, we calculate the damping length averaged along the 
spicule, Ld, as 



- _ \ r" 



(66) 



where s represents the direction along the spicule and H is the 
height of the chromosphere above the photosphere. We take 
H - 3,000 km. Equation (l66T l is used to calculate the averaged 
values of both Ld.ra and Ld.in- We plot in Figure |4fb) the av- 
eraged values of the damping lengths as functions of the wave 
period. First, we obtain that the averaged Ld.in is several or- 
ders of magnitude longer than the averaged Ld.ra in the range 
of periods taken into account in Figure EJb). This means that 
ion-neutral collisions have little impact on wave propagation. 
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On the contrary, the averaged Ld.ra for periods less than 10 s 
is smaller than or of the same order as the height of the chromo- 
sphere. This result points out that only waves with periods longer 
than 10 s are able to reach the coronal level in the form of kink 
motions. Waves with shorter periods reach the corona as small- 
scale azimuthal motions, which are unobservable with present 
day instruments. This effectively imposes a lower limit for the 
period of kink waves observable in the corona. T his is consistent 
with the observed periods of cororial waves (e.g^lTomczyk et all 
l2007HTomczvk & McIntoshll2009HMcIntosh et al. 11201 Hk 

6. Conclusion 

In this paper we have investigated resonant Alfven waves in par- 
tially ionized plasmas. We find that the conserved quantity at 
the resonance and the jump of the perturbations across the res- 
onant layer are the same as in fully ionized, ideal plasmas. We 
have derived expressions for the damping lengths due to reso- 
nant absorption and due to ion-neutral collisions for the case of 
propagating kink waves in straight magnetic tubes. In the limit of 
large collision frequencies, the damping length due to resonant 
absorption is inversell y proportional to the frequency as in the 
fully ionized case (seel Terradas et al]|20lQi) . whereas the damp- 
ing length due to ion-neutral colhsions is inverselly proportional 
to the square of the frequency. We have applied the theory to the 
case of chromospheric kink waves propagating from the photo- 
sphere to the corona. We conclude that the solar chromosphere 
acts as a filter for kink waves. Waves with peridos shorter than 
10 s, approximetely, reach the corona in the form of small-scale 
azimuthal waves. Only waves with periods longer than 10 s can 
be observed in the corona as kink waves. 
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